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ABSTRACT 

We present detections of the ground-state lio loi transition of ortho- 
H2O at 557 GHz in 18 molecular outflows based on data from the Submillimeter 
Wave Astronomy Satellite (SWAS). These results are combined with ground- 
based observations of the J=l-0 transitions of ^^CO and ^^CO obtained at the 
Five College Radio Astronomy Observatory (FCRAO). Data from Infrared Space 
Observatory (ISO) for a subset of the outflows are also discussed. Assuming the 
SWAS water line emission originates from the same gas traced by CO emission, 
we find that the outfiowing gas in most outfiows has an ortho-H20 abundance 
relative to H2 of between about 10~^ and 10~®. Analysis of the water abundance 
as a function of outflow velocity reveals a strong dependence. The abundance 
of ortho-H20 increases with velocity and at the highest outflow velocities some 
of the outflows have relative ortho-H20 abundances of order 10~^. However the 
mass of very high velocity gas with such elevated H2O abundances represents 
less that 1% of the total outflow gas mass. The ISO LWS observations of high-J 
rotational lines of CO and the 179.5 //m transition of ortho-H20 provide evidence 
for a warmer outflow component than required to produce either the SWAS or 
FCRAO lines. The ISO line flux ratios can be reproduced with C-shock models 
with shock velocities of order 25 km s"^ and preshock densities of order 10^ cm~^; 
these C-shocks have post-shock relative water abundances greater than 10~^. The 
mass associated with the ISO emission is also quite small compared with the total 
outflow mass, and is similar to that responsible for the highest velocity water 
emission detected by SWAS. Although the gas responsible for the ISO emission 
has elevated levels of water, the bulk of the outflowing gas has an abundance 
of ortho-H20 well below what would be expected if the gas has passed through 
a C-shock with shock velocities greater than 10 km s^^. Gas-phase water can 
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be depleted in the post-shock gas due to freeze-out onto grain mantles, however 
the rate of freeze-out is too slow to explain our results. Therefore we believe 
that only a small fraction of the outflowing molecular gas has passed through 
shocks strong enough to fully convert the gas-phase oxygen to water. This result 
has implications for the acceleration mechanism of the molecular gas in these 
outflows. 

Subject headings: ISM: jets and outflows — ISM: molecules — ISM: abundances — 
Stars: formation — Stars: winds, outflows 
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Introduction 



Stellar winds are believed to play a pivotal role in the process of star formation. During 
the accretion phase of star formation, stellar winds carry away excess angular momentum 
that would otherwise prevent further collapse. The winds interact with and accelerate the 
surrounding medium, and as a result, drive molecular outflows with masses often much 
greater than that of the young star itself. However, the exact mechanism by which these 
winds interact with and accelerate the surrounding medium is still a matter of debate. 
A number of different models have been proposed including iets with bow shocks, jets 



with turbulent entrainment, and wide angle winds (lArce fc Goodman 



2OO2I). The various 



models result in different spatial and velocity distributions of the entrained material, and 



therefore one approach to distinguish between the pro posed mode 



s is t o map the spatial 



20061 ). Alternatively, the 



and velocity structure of the outflowing molecular gas (lArce et al.l 
various acceleration mechanisms heat the gas in different ways, resulting in changes in the 
chemistry of the outflow as ices are sublimated and endothermic reactions pathways are 
opened. Thus, determinations of the chemical abundances of the outflowing m olecular gas 



may p ermit us to distinguish between the various mechanisms. The review by 



ArceetaL 



( I2OO6I ) suggested that molecules such as SiO, CH3OH, H2O and s ulfur-bearing rnolecule s 



Bachiller et al. 



(1200lh 



have their abundances affected by outflow activity. For example, 
study of LI 157 revealed significant abundance variations and suggested that the chemistry 
of the outflows may be useful in understanding the time evolution of outflows. 



The abundan ce of water can be strongly affected by shocks (IKaufman &: Neufeld 



1996 



Bergin et al. 



19981 ) such as those predicted to exist in outflow regions. Moderate velocity 



shocks are capable of producing large abundances of water by releasing frozen water from 



dust grain s and by driving all 



reactions ( Kaufman fc Neufeld 



ree ox ygen into water through a series of gas phase chemical 



19961 ). Thus, water should be a good tracer of shocks. 
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The chemical reactions that produce water have been thoroughly modeled and it has 
been shown t hat the production mechanism and resulting abu ndance is very sensitive to 



temperature (lElitzur &: de Jong 



1978 



Elitzur fc Watson 



19781 ). At gas temperatures below 



~300 K water is formed most readily through a series of ion molecule reactions that lead 
to the formation of HsO"^. The dissociative recombination of HsO^ has several possible 
outcomes, b ut has been measur ed in the lab to produce water in approximately 25% of the 



interactions (jJensen et al. 



20001 ) ■ This process of water production is relatively slow and 
modelling has shown that in a quiescent medium of density 10^ cm~'^ and temperature 30 



K, a water abundance o 



years (IBergin et al. 



19981). 



roughly 10 ^ relative to molecular hydrogen is achieved after 10^ 



Once the temperature of the gas rises above ~300 K a series of endothermic neutral- 
neutral reactions are activated. At elevated temperatures, these reactions are very rapid 
compared to the ion mo l ecule reactions, arid all fr ee oxygen is quickly d riven into water 



( lElitzur &: de Jong 



1978 



Elitzur &: Watson 



19781) 



Kaufman fc Neufeldl ( 119961 ) show that 



the passage of a moderate C-shock with a velocity greater than ~10 km s ^ is able to elevate 
temperatures above 300 K long enough to produce gr eatly enhanced wate r abundances on 



Bergin et al. 



1998|) showed that the 



the order of ~ 10~^ relative to H2. More importantly, 
enhanced water abundance persists in the post-shock gas after the gas has cooled and to 
significantly reduce the water abundance required ~ 10^ years even at densities as high as 
10^ cm~'^. Since outflow ages are typically estimated to be of order 10^ years or less, any of 
the outflowing molecular gas that has passed through a moderate shock should be imprinted 
with this greatly elevated water abundance. After ~ 10^ years the abundance of water 
will drop as the water chemistry comes into equilibrium at the lower gas temperature and 
water freezes onto the cold dust grains. Evidence for this enhanceme nt is quite apparent 



in the water qbseryations of some outflo ws such as in Orion-BN/KL (iHarwit et al 



Wright et al 



200G 



Melnick et al 



2000bh 



1998 : 
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The ideal nature of water as a tracer of shocked regions is countered by the difficulty 
of detection from ground based observatories. This problem was overcome by the successful 
launch of three satellites, the Infrared Space Observatory (ISO), the Suhmillimeter Wave 
Astronomy Satellite (SWAS) and Odin. ISO was cap able of detecting a number of water 



transitions excited in gas warme r than about 80 K 



flMelnick et al. 



2000al ) and Odm (iNordh et all 



2003 



fjClegg et al 



Hjalmarson et al 



199fih wh ile SWAS 



20031 ) can observe 



the fundamental ortho-H20 lio — loi transition at 538.3 /im. More recently, the Spitzer 
Space Telescope has detected several mid-infrared water transitions toward NGC 2071 



(IMelnick et al. 



20071 ). The fundamental transition of ortho-H20 can be readily excited at 
temperatures greater than about 20 K, permitting SWAS and Odin to trace out much cooler 
post-shock gas than ISO or Spitzer. SWAS has detected ortho-H20 emission associated 
with the quiescent dense molecular gas in numerous cloud cores. The relative ab undance of 
ortho- H^O to H2 in these regions was found to be on the order of ~10^^ to 10^^ (ISnell et al. 



2000bl ). roughly 2 orders of magnitude less abundant than predicted by quiescient gas phase 



chemist ry. The discrepenc y is likely the result of water ice mantle formation on the dust 



grains (IBergin et al. 



2000). 



In this paper, we present SWAS detections of H2O emission from 18 well studied 
outflow regions. We selected nearby outflows (within 1 Kpc) with significant SWAS 
observations that we were able to obtain additional mapping data from the Five College 
Radio Astronomy Observatory (FCRAO). These outfiows are driven by a mix of high and 
low luminosity Young Stellar Objects (YSOs). Eight of these outfiows (L1448-mm, NGC 
1333-SVS13, HH25mm, p Oph A (VLA1623), L1689N, Ser SMMl, L1157, L1228) are driven 
by relatively low luminosity YSOs that are in most cases individual low mass stars. The 
other ten sources (GL490, Orion KL, OMC 2, NGC2071, MonR2, NGC2264 D, NGC2264 
C, IC 1396N, S140, Ceph A HW2) are driven by much higher luminosity YSOs that are 
associated with newly forming groups and clusters. We combine the SWAS II2O detections 
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with maps of ^^CO and ^^CO emission obtained at the FCRAO 14-m telescope to determine 
the water abundance in the outflowing gas. Additionally, we use published ISO results for 
six of these outflows to further expand our understanding of the SWAS results. 



2. Observations 



2.1. SWAS 



SW AS was a NASA Small Explorer Mission that operated successfully from 1998 
to 2005 (IToUs et al.l 12004 ) . SWAS simultaneously observed the lio — loi transition of 
ortho-HsO at 556.936 GHz, the J = 5-4 transition of ^^CO at 550.926 GHz, the 3,1-3,2 
transition of O2 at 487.249 GHz, and the ^Pi Pq transition of [CI] at 492.161 GHz. 
O2 emission was not detected toward these outflows and these observations will not be 
discussed further. The high spectral resolution of SWAS (~0.6 km s~^) permits the 
kinematic separation of the outflowing gas from the ambient material within star forming 
regions. However, a detailed analysis of the distribution of emission in these outflows is 
hindered by the limited angular resolution afforded by the small aperture of SWAS. The 
SWAS beam is elliptical, and at the frequency of the H2O and ^^CO transitions has angular 
dimensions of 3'. 3 x 4'. 5 and at the frequency of the O2 and [CI] transitions has angular 



dimensions of 3'. 5 x 5'.0 (IMelnick et al. 



2000al ). Only one pointing was obtained for each 



source; however, in most outflows the large beam encloses the entire region of outflow 
activity. 

Observations of these outflows were obtained by SWAS over the entire period of mission 
operations. The positions observed for the 18 outflows are listed in Table 1. The data was 
acquired by nodding the satellite alternately between the source position and a reference 
position relatively free of ^^CO emission. Each position was observed over many orbits and 
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the total integration time on source is listed in Table 1. For each observation, an equal time 
was spent integrating off source toward the reference position. The on-source integration 
time varied from about 6 to 90 hours and thus, the rms noise level achieved varied from 
source to source. The shortest integration time was for 0MC2, where the rms noise 
measured in the baseline for the H2O and ^^CO spectra was 0.028 K. The noise decreased 
as the square root of time and thus for one of the longer integration sources, p Oph A, an 
rms of 0.008 K was obtained. For the source with the longest integration time, GL490, the 
baselines were poorer and the rms noise obtained was only 0.010 K. Details concerning dat a 
acq uisition, calibratio n, and reduction with SWAS are presented in 



Melnick et al. 



and 



Tolls et al. 



(I2000af ) 



(120041 ). The data sho wn in this paper are not corrected for the measured 



SWAS main beam efficiency of 0.90 (IMelnick et al. 
been applied in our analysis. 



2000al ): however, this correction has 



2.2. FCRAO 

During the fall of 2003 and the spring of 2004, the FCRAO 14-m telescope was used to 
obtain maps of the ^^CO emission at a rest frequency of 115.271 GHz and the ^^CO emission 
at a rest frequency of 110.201 GHz. The map centers were the same positions as those used 
for the SWAS observ ations. For all observations, the 32-pixel SEQUOIA array receiver 



( lErickson et al. 



19991 ) was used and data was obtained using an On-The-Fly observing 
technique. For each source the observations were resampled to form maps approximately 
6'.0 X 6'.0 in Aa and A(5 in extent with data spaced by 25". The spectrometer for each 
pixel was a digital autocorrelator with a bandwidth of 50 MHz and 1024 spectral channels 
per pixel leading to a channel spacing of approximately 0.13 km s~^ for both ^^CO and 
^^CO. The full-width-at-half-maximum beam size of the FCRAO telescope at the ^^CO 
frequency is 45". The spectra in these maps typically have an rms noise measured in the 
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baseline of 0.06 K for ^^CO and 0.17 K for ^^CO. The two exceptions are L1689N and p 
Oph A which transit at low elevation at the latitude of the FCRAO site and had an rms 
noise measured in the baseline of only 0.14 K for ^^CO and 0.47 K for ^^CO. The CO data 
shown in this paper have not been corrected for main beam efficiency, estimated to be 
0.45 at the frequency of the ^^CO line and 0.49 at the frequency of the ^'^CO line. These 
corrections have been applied in all subsequent analysis. 



2.3. ISO 



ISO had a much greater wavelength coverage than SWAS and was capable of measuring 
the far-infrared emission from a va riety of ortho- an d para-H20 transitions as well as a 



series of high-J CO emission lines ( IClegg et al. 



19961 ). Pubhshed ISO res ults are aval 



for seyen o f the outflows in our sample: these outflo ws are L 1448- mm (jNisini et al. 



1999 



200G 



HW2 ( 



Froebrich et al 



Froebrich et al 



(ICeccarelli et al 



and Ser SMMl (iLarsson et al. 



2002), N GC1333 SVS13 flMohnari et al 



2002), LI 15 7 flNeufeld et al 



1998). HH25mm (Benedettini et al. 



20001: 



Giannini et al 



able 



2000h. C eph A 



200lh. L16 89N 



2000), Orion KL flLerate et al.ll2006h . 



20021 ). Water lines were detected by ISO in all of these 
outflows except Ceph A HW2. The published line fluxes that we use in this paper were 
all obtained with the LWS instrument. With the exception of Orion KL, these spectra 
were obtained at relatively low spectral resolution (R ~ 300) so they contain little velocity 
information. The lowest lying water line observed by ISO is the 2i2 — loi transition at 179.5 
yUm. This transition originates from a higher energy level (E/k = 80 K) than the transition 
observed by SWAS and therefore may probe a higher temperature component of water. 
ISO has a smaller beam than SWAS (~75") but the resolution is still too poor to permit a 
detailed spatial analysis. 
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3. Summary of Observations 



In Figure 1, we present tlie SWAS spectra of tlie lio — loi transition of ortho-H20 
for tlie 18 outflows in our sample. Each spectrum has only had a linear baseline removed 
Seven of the outflows p resented in this pap er have had SWAS data pu blished previously : 



NGC 2071 and L1157 flNeufeld et al. 



20021), NGC1 333 SVS13 



and L1689N fIStark et al. 



( Bergin et al 



20001). L1448-mm and IC13 96N flBenedettini et al. 



20031 ). S140 and Mon R2 (IBoonman et al. 



20031) 



20M). L1157, 1C1396N and S140 have been reobserved by SWAS 



since the time of these publications and the complete spectra are shown in Figure 1. The 
initial observations of GL490 had a ripple in the baseline and so this source was reobserved. 
For the new observations the local oscillator setting was shifted to a value that limited the 
extent of useable spectrum on the blue side of the line. With the hmited spectral coverage, 
placement of the baseline in GL490 is more questionable than in the other outflows, and 
therefore the extent of the outflow emission in the blue wing is very uncertain. Finally, for 
L 1448- mm the position observed by SWAS and the center of the CO mapping at FCRAO 
was offset ~3 ' north-west of L1448-mm, the center of the large bipolar outflow. Thus, our 
observations of L 1448- mm primarily probe the blueshifted lobe of the outflow. The SWAS 
results for ^^CO and [CI] will be presented later in the paper. 

We have obtained nearly complete CO maps of all 18 outflows. In only L1228, 
L1448-mm and NGC2264 D does the molecular outflow extend appreciably beyond the 
6' X 6' region mapped. The CO maps of 111125mm and Ser SMMl show obv ious evidence of 



multiple outflows in the region mapped, and this has been noted previously (iGibb &: Heaton 



1993 



Davis et al. 



19991 ). The SWAS spectra, although centered on HH25mm and Ser 
SMMl, include emission from other outflows in these regions. Since all 18 of these outflow 
regions have been extensively studied and there are many published papers with outflow 
maps, we do not believe our CO maps provide sufficient new information to justify showing 
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here. A summary of the pro perties of t 
in the catalog introduced in 



WuetaL 



lese outflows and their driving sources are provided 



(120041 ) 



We convolved the CO data to the SWAS resolution to compare H2O and CO line 
profiles. However, since the outflows do not fill the SWAS beam, convolving the CO data to 
match the large SWAS beam substantially degraded the signal to noise in the high velocity 
CO line wings making profile comparison difficult. So instead, we have chosen to co-average 
only those ^^CO spectra which have an integrated intensity in either the red or blue outflow 
wings that is at least one-half of the integrated intensity of the peak wing emission. In 
these same positions we have also co-averaged the ^^CO spectra. Both ^^CO and ^^CO 
spectra averaged in this manner are shown accompanying the water spectra in Figure 1. It 
is worth noting a few features in the averaged CO spectra. First, in Ceph A IIW2 there 
is a additional velocity feature at Visr = 5 km s^^ that is unlikely related to the outflow. 
In GL490 and L1228 a weak negative feature appears in these spectra, suggesting that the 
reference positions were not free of ^^CO emission at all velocities. We excluded the velocity 
intervals associated with these features in the three sources in our subsequent analysis. 

A cursory inspection of the CO and H2O line profiles reveal many similarities that 
suggest that the emission in these lines likely originates in the same outflow component. 
With the exceptions of Mon R2, NGC 2264 D, and S140, the H2O emission can be traced 
to a much higher outflow velocities than the CO emission. This effect is most pronounced 
in the low luminosity outflows, such as LI 157 or 111125mm, and may be due to differences 
in the optical depth of these lines. The CO J=l-0 transition is intrinsically weak, and 
we might expect the optical depth of this line to be less than the water line. Higher 
rotational transitions of CO, which have larger optical depths, do reveal higher velocity 
emission. A good example is L1448-mm where obervations of the J=2-l transition of CO 



( iBachiller fc Cernicharo 



1990) reveal a velocity extent more similar to what is seen in the 
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water line. Although the detection of outflows is more difficult in the J=l-0 transition, 
the smaller optical depths make it ideal for deriving the outflow column density. We also 
find that the intensity ratio of the ortho-H20 emission to ^^CO emission increases with 
increasing outflow velocity. Later in the paper we will model these emissions, including the 
effects of optical depth, and test whether the observed changes in the intensity ratio is due 
to variations in the water abundance with velocity. 

Finally we note that nearly all of the H2O spectra show evidence for self-absorption 
similar to that seen in ^^CO. The self- absorption feature is often much less evident in the 
co-added ^^CO spectra shown in Figure 1 than it is in the individual ^^CO spectra obtained 
toward the center of the outflow. The self-absorption feature seen in both H2O and ^^CO 
presumably arises due to foreground, quiescent, low-excitation gas absorbing the central 
part of the broad emission produced by the outflowing gas. 



4. Analysis and Results 

To derive the water abundance from the SWAS ortho-H20 observations, we must 
quantify the density, temperature and column density of the outflowing molecular gas. 
Densities and temperatures are extremely difficult to determine in molecular outflows 



and few studies have attempted to measured these properties. Some results ar e aval 
concerning the gas temperature in outflows and are summarized in 



WuetaL 



(12004 ) 



able 



However for the density, it is difficult to find any reliable estimates. In this paper we 
will simply assume values of temperature and density typical of star forming cores. 
Specifically, we set the density to 10^ cm~'^ and the temperature we assume varies from 
30-100 K, depending on the luminosity of the driving source. For outflows associated 
with low-luminosity YSOs (L1448-mm, NGC 1333 SVS13, HH25mm, p Oph A, L1689N, 
Ser SMMl, L1157, L1228 and IC1396N) we assume a temperature of 30 K. For outflows 



- 13 - 



associated with more luminous YSOs ( GL490, 0MC2, NGC2071, MonR2, NGC 2264 D, 
NGC2264 C, S140 and Ceph A HW2 ) we assume a temperature of 50 K and for Orion 
KL we assume a temperature of 100 K. These assumptions are an oversimphfication of the 



complex temperature and density structure like^ 
of outflows were systematically investigated by 



y to be pr e sent i n these outflows. A number 



Levreault 



( 1l988l ) and the range of densities 



and temperatures estimated for outflows in this study are consistent with our assumptions. 

Higher temperature gas is present in nearly all of the outflows studied here as evidenced 
by the presence of shock excited optical emission (Herbig-Haro objects) as well as highly 
excited lines of a variety of molecular species seen at infrared and far infrared wavelengths. 
In particular, ISO has detected high-rotational transitions of CO and H2O. However, the 
bulk of the outflowing molecular gas likely has much lower gas temperatures. In our analysis 
we will first address the emission from the cooler gas (T = 30-100 K) that we believe 
predominately gives rise to the emission in the lowest rotational transitions of CO. In § 4.6 
we will address specifically the emission from the warmer gas within the outflows. 



4.1. Outflow Gas Column Density 

We derive the column density of outflowing gas based on the ^^CO J=l-0 emission 
assuming the outflow gas temperature and density as described in the previous section. 
It is necessary to begin by identifying the velocity intervals over which emission from 
the outflowing gas dominates the ambient cloud emission. As our guide for setting these 
intervals, we have used the observed ^^CO spectra and defined the break in line profile 
shape as the division between ambient and outflow emission. The intervals so defined are 
given in Table 2 and illustrated in Figure 1. Based on these velocity intervals, we compute 
the integrated intensity of the CO emission in each map position for each outflow. The 
emission from the outflowing gas within the velocity interval around the quiescent cloud 
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velocity is excluded in our analysis; however, this velocity interval is relatively small and its 
exclusion will not strongly affect the results. 

The integrated intensity of CO can be converted to a column density, assuming that 
the emission is optically thin and in LTE. At the high densities assumed for the molecular 
outflow gas (much greater than the critical density for ^^CO) the assumption of LTE is 
extremely good. The total CO column density can be written as 



N = 1.15 X 10^^e^-^^/^(^ + 0.36r) j T^bdv, 



where T is the excitation temperature (which is equivalent to the gas temperature in LTE), 
and Tmb is the main beam antenna temperature (antenna temperature corrected by the 
main beam efficiency) of the J=l-0 transition of CO. 

Although our assumption of LTE is likely valid, we have no reason to believe that the 
observed high velocity ^^CO emission is optically thin. Our simultaneous observations of the 
^^CO emission from the outflows enable us to test this thin assumption and, if necessary, 
correct for any optical depth effect. We estimate the optical depth of this gas by comparing 
the ob served ratio of ^^CO/^^C O emission to the assumed local interstellar medium ratio 



1990l ). Assuming equal excitation temperatures for ^^CO and 



of 65 (iLanger fc PenziasI 
^^CO, the relation between the observed isotopic line ratio, R, and the optical depth can be 



written as 



1 - e--/65 ' 



where r is the optical depth of the ^^CO emission. If we further assume that the ^^CO 
emission in the high velocity gas is optically thin (r/65 <C 1) then the ratio becomes 

i?^^(l_e-). 

r 
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We then apply a correction to the column densities derived assuming the emission is 
optically thin, and the corrected column density is then simply 



The sensitivity of our ^^CO maps is insufficient to detect high velocity gas in individual 
spectra. Because of the limited angular extent of the outflows, if we average over the entire 
map or even over the SWAS beam, we severely reduce the outflow signal in the averaged 
spectrum. Therefore, we have averaged only those spectra that show evidence of the outflow 
in ^^CO emission, following the method used to construct the ^^CO and ^^CO spectra 
shown in Figure 1. The redshifted and blueshifted gas were considered separately, unlike 
the method used to produce the spectra shown in Figure 1, and only those positions with 
^^CO integrated areas greater than one-half the peak value detected for the redshifted or 
blueshifted outflow in that source were selected and averaged together. The ^^CO profiles 
were created by averaging the same positions selected by the ^^CO emission. The observed 
isotope ratio in either the redshifted or blueshifted outfiow is then measured from these 
averaged spectra. The measured isotopic ratios are presented in Table 2. The measured 
value of R is typically between 20 and 30 which correspond to optical depths of 3.0 to 
1.8 and to column density corrections of a factor of 3.3 to 2.2. The average optical depth 
correction is used to correct the optically thin derived column density throughout each 
outfiow lobe. 



Nthick — Nthin 



T 




4.2. 



CO J=5-4 Emission - Test of Outfiow Model 



The SWAS observations of the J=5-4 transition of ^^CO can be used to test our adopted 
phyical properties and derived gas column density of these outfiows. The emission in this 
line originates from an upper state with an energy of E/k = 79 K; thus it is particularly 
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sensitive to the assumed temperature in the range considered for these outflows. Signiflcant 
i^co J=5-4 emission was detected in seven of the outflows (Ceph A HW2, GL490, MonR2, 
NGC2071, Orion KL, OMC 2 and S140) aU driven by luminous YSOs. For outflows 
associated with lower luminosity YSOs, ^^CO J=5-4 emission was either extremely weak or 
not detected. This general trend is consistent with our assumption that the outflowing gas 
associated with more luminous YSOs is warmer than that in the lower luminosity YSOs. 
However, a more quantitative test of our physical model can be made. 

Using the computed CO gas column density and the assumed temperature and density 
for the outflows, we can predict the emission in the ^^CO J— 5-4 transition at each location 
in the outflow where we have determined the outflow column density. As in our previous 
column density determination, we assume that the abundance ratio of CO to ^^CO is 
65. We convolve the map of predicted emission with the SWAS beam, and compute the 
predicted ^^CO J=5-4 hue flux in each of the outflows and compare that with the line 
flux observed by SWAS. We flnd that for most of the sources with signiflcant ^^CO J=5-4 
emission that the modeled and observed line fluxes are consistent within a factor of 2. The 
agreement is remarkably good, since only small changes in the density and temperature can 
greatly affect the hue flux. The most discrepant outflows were Mon R2 and GL490, where 
the modeled line flux was 3 to 6 times larger than that observed. In three of the outflows 
associated with luminous YSOs, IC1396N, NGC2264 C and NGC2264 D, there was no 
detectable ^^CO J=5-4 emission; however, in all cases the non-detections are consistent 
with the model predictions. For the outflows associated with low luminosity YSOs, only 
p Oph A and L1689N have possible weak outflow detections in the ^^CO J=5-4 line, and 
the other sources (L1448-mm, NGC1333-SVS13 SVS3, HH25mm, Ser SMMl, L1157, and 
L1228) were not detected. We flnd again that for these low-luminosity outflows with weak 
detections or only upper limits the results are consistent in most cases with our model. The 
one exception is HH25mm, where the modeled flux exceeds the 3a upper limit by a factor 
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of approximately 3 to 5. In all cases where there is poor agreement, the model predictions 
exceed the observed line flux suggesting that we have either overestimated the temperature 
or density in our outflow model. 

We have used two of the most discrepant outflows, MonR2 and HH25mm, to estimate 
the magnitude to which we may have overestimated cither the density or temperature 
in our outflow models. Unfortunately with only the ^^CO J=5-4 data it is impossible to 
derive a unique temperature and density for the outflows, since both parameters affect 
the observed line flux. For these calculations we have recomputed the gas column density 
using the modifled temperature or density and then recomputed the predicted ^^CO J— 5-4 
emission. In MonR2, we have found that by reducing the temperature from 50 to 30 K we 
can obtain good line flux agreement. However, we can also maintain the temperature at 50 
K and reduce the density by a factor of 5 and obtain good agreement. Of course the correct 
solution could be a combination of both a lower temperature and a lower density. Likewise 
in HH25mm, we can achieve good flux agreement by either lowering the temperature from 
30 to 20 K or by decreasing the density by a factor of 4. Nevertheless, in most outflows 
the line flux agreement suggests that the temperature is probably accurate to 20% and the 
density accurate to within a factor of two. 

4.3. Mass of Outflows 

The mass of outflowing redshifted and blueshifted gas in each outflow can be found 
using our map of the CO column density and the distances given in Table 1. We have 
assumed a CO/H2 ratio of lxlO~^ for the mass determinations which are summarized in 
Table 2. The statistical errors on the observed integrated intensity of ^^CO are small, and 
the uncertanity in our mass estimates are dominated by the uncertainties in our assumed 
temperature, correction for optical depth, and CO to H2 conversion. Uncertainties in 
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density have little impact on the mass determination, unless it is substantially lower that 
the value we assumed. 

Our method for estimating the water abundance requires a knowledge of the spatial 
distribution of the column density of outflowing gas. For this analysis, which we describe 
in the next section, we have only included positions which have a 3cr detection of ^^CO. 
The mass estimates in Table 2 are based only on positions within the outflow with CO 
detections, and thus may underestimate the total outflow mass. Statistically we can provide 
a much better estimate of the total outflow mass by using all the positions, regardless of 
whether the integrated intensity is positive or negative. We have recomputed the mass of 
each outflow without applying a detection threshold, and find masses in the redshifted and 
blueshifted outflows that are at most only ten percent larger than those presented in Table 
2. 



4.4. H2O Abundance 

We first estimate the ortho-H20 abundance in the outflows assuming that the H2O 
emission seen by SWAS originates in the same gas as traced by the CO J=l-0 emission. 
Later we will explore whether this is reasonable assumption and if there are other 
possibilities for its origin. The integrated intensity of the SWAS H2O line in each outflow 
has been determined using the same velocity intervals that were used for CO. The observed 
integrated intensity and line flux for the redshifted and blueshifted emission is given in 
Table 2. We model the water emission using a statistical equilibrium code that uses the 
Large Velocity Gradient (LVG) approximation to account for radiation trapping. The LVG 
assumption should be a valid approximation for the obser ved broad-line emissi on from 



these outflows. Collisional rate coefficients are taken from 



Phillips et al 



19961 ). Since 



the ortho-H2 and para-H2 collisional rate coefficients with H2O are different by nearly 
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an order of magnitude, the assumed ratio of ortho-H2 to para-H2 is important. We have 
assumed that t he ra tio of ortho- to para-H2 is in LTE as found for the outflows studied by 



Neufeld et al 



(120061 ) . Thus for our assumed gas temperatures of 30, 50 and 100 K, the 
ratio of ortho- to para-H2 is approximately 0.03, 0.3 and 1.6, respectively. We include the 
five lowest levels of ortho-H20 in our calculation. 

The physical inputs to our H2O model are the temperature, density, and column 
density distribution of the outflowing gas. We have used the same temperature and density 
as described earlier for computing the CO column density. The process of solving for the 
abundance begins by assuming a water abundance and then computing the predicted water 
emission at each position in the outflow. We convolve the predicted emission with the 
SWAS beam, and then compare with the observed integrated intensity. The abundance 
of ortho-H20 is then varied until we achieve agreement between model and observations. 



This approach is very similar to that used by 



Snell et al 



(l2000al ) and has the advantage of 



correctly accounting for optical depth variations and beam dilution effects across the large 
SWAS beam based on our measured gas column density distribution. Since the input to our 
model is the measured CO column density, our analysis provides a direct determination of 
the ortho-H20 abundance relative to CO. We infer an abundance of ortho-H20 relative to 
H2 by assuming a ^^CO to H2 ratio of 10~^. The resulting ortho-H20 abundances relative 
to H2 are presented in Table 2. 

We have investigated the impact on the relative ortho-H20 abundance for different 
assumed temperatures and densities. Unless the assumed density is dramatically decreased, 
density uncertainties have little impact on the determination of the CO column density. 
However, density does have a big effect on ortho-H20 abundanc e. Since the water e mission 



in these outflows is approximately in the effectively thin limit ( iNeufeld et al. 



2000), the 



abundance derived will be inversely proportional to the assumed density. Thus, if we were 
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to assume a density a factor of two larger, the resulting ortlio-H20 abundance will be a 
factor of two smaller. Changes in the assumed temperature, on the other hand, will effect 
both the derived column density of CO and ortho-H20 abundance. We have used the 
redshifted outflows of MonR2 and L1689N as test cases. In each case we have increased 
and decreased the temperature by a factor of two and recomputed the CO column density 
and ortho-H20 abundance. We found for these outflows that increasing the temperature by 
a factor of two results in a decrease in the relative abundance of ortho-H20 by a factor of 
3 in Mon R2 and by a factor of 8 in L1689N. Decreasing the temperature by a factor of 
two results in an increase in the derived ortho-H20 abundance by a factor of 9 in Mon R2 
and by a factor of 5 in L1689N. Thus, our derived ortho-H20 abundances are very sensitive 
to the assumed physical conditions in the outflowing gas, which is unfortunate since these 
conditions are very poorly known. Thus, the systematic errors introduced by our poor 
knowledge of the outflow density and temperature dominate over the uncertainties in the 
H2O line fluxes in determinations of the water abundance. 

The abundance of H2O relative to H2 spans a wide range of values, from 10^^ to 10^®. 
This range of abundances is so large, we believe it is unlikely due to only uncertainites 
in the physical conditions and must reflect true abundance variations. We note that the 
results for NGC2071 and L11 57 are in good agreement with the abundances reported 



earlier by 



Neufeld et al 



(I2OOOI ) that are also based on the SWAS data. We find the highest 



H2O abundances are found in Orion KL and LI 157, and the lowest abundance in Mon R2. 



Our derived ortho-H20 abundances in nearly all o: 
derived for quiescent molecular cloud gas (~10^^; 



the outflows are much larger than that 



Snell et al. 



2000bl ) which suggests some 
enhancement of the water abundance in outflows. The largest abundance derived, ~10~^ 
in Orion KL, is still less than one might expect if all of the free oxygen had been coverted 
into water as is expected for a moderate shock. It is very important to remember that 
these abundance estimates are averages over the entire SWAS beam, and do not rule out 
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the possibility of much higher abundances existing over a small portion of these outflows. 
In fact, in the warm gas probed by ISO, the abundance of H^Q i s as high as 5x10^ "^ 



1998 : 



Giannini et al. 



200l|). 



in the warm gas in Orion KL and L1448-mm (iHarwit et al. 
The relationship of the warm gas probed by ISO and the emission seen by SWAS will be 
examined in more detail in § 4.6. 



4.5. Velocity Dependance of H2O Abundance 

The high spectral resolution of SWAS enables us to further investigate the abundance 
of H2O as a function of velocity. In the previous section we assumed that the ortho-Il20 
abundance is constant thoughout the outflow. We will now investigate whether the 
abundance of ortho-H20 relative to ^^CO varies with outflow velocity. 

We have divided the redshifted and blueshifted emission from the outflowing gas in 
CO and H2O into velocity bins, each of width 5 km s^^, continuing to higher velocities 
until ^^CO is not detected. For a velocity interval to be analyzed, at least two positions 
within the outflow must have significant detections of ^^CO emission. The integrated H2O 
intensity and line flux as a function of velocity are presented in Table 3 for each outflow. 
Optical depth corrections were considered separately for each velocity interval using the 
same averaged CO spectra as described before. Any interval with a non-detection of 
^^CO was assumed to be optically thin, which is often the case at the higher velocities in 
the outflows. We then used the same technique as discussed above to model each of the 
individual velocity intervals. The results of our analysis are presented in Table 3, which 
includes the observed CO isotopic ratio, the mass of gas, and the ortho-Il20 abundance as 
a function of velocity. 

Our abundance analysis requires ^^CO detections. Although, earlier we suggested that 
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including only 3 a detections had only small effect on the total outflow mass, the impact 
of truncation is more severe in the higher velocity bins where the signal to noise of the 
^^CO emission is relatively low. We have recomputed the mass in each velocity bin without 
applying a detection threshold, and found that in the most discrepant bins, which are at 
the highest velocities, the mass could be underestimated by a factor of 2 to 3. Thus, the 
ortho-H20 abundances quoted may by overestimated by a similar factor. Despite these 
uncertainties, we find that in every outflow the derived abundance of ortho-H20 increases 
steadily with increasing outflow velocity. Even outflows that did not show obvious disparate 
velocity extents between CO and H2O show ortho-H20 abundance variations with velocity. 
Some of the most pronouced abundance variations with velocity are in the outflows driven 
by low-luminosity YSOs; for example in NGC1333-SVS13 the ortho-H20 abundance in both 
the redshifted and blueshifted high velocity gas increases by two orders of magnitude from 
the lowest velocity outflowing gas to the highest velocity outflowing gas. At the highest 
outflow velocities, the abundance of ortho-H20 relative to H2 gas can be as high as 1 x 10~^, 
approaching values that one might expect if all of the oxygen not in CO is converted 
into H2O as anticipated if the gas has recently passed through a moderate velocity shock. 
However, since most of the outflow mass is at relatively low outflow velocities, the mass 
of gas with this highly elevated ortho-H20 abundance represents less than 1% of the total 
mass of outflowing gas, although, it accounts for a much larger fraction of total water line 
flux. 

4.6. ISO Analysis 



ISO with the LWS detected both high-J CO lines and a number of ortho- and para-H20 
hues in eight of the molecular outflows in our survey. The lines observed by ISO likely 
arise in a much warmer gas than is responsible for the bulk of the low-J CO outflow 
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emission that we have modeled. For water, the lowest excitation line observed by ISO (the 
2i2 — loi transition at 179.5 /xm) has a flux in most sources more than 10 times larger 
than that predicted by the outflow models discussed in previous sections. Thus, the ISO 



water emission must be produ ced by a warmer, 



outflowing gas 



Molinari et al. 



Many authors fjHarwit 



2000 



et al 



Benedettini et al.l 



20001 : 



1998 



i kely shock-excited component, of 



the 



Ceccare 



Giannini et al. 



i et al. 



1998 



Nisini et al. 



20001 : 



20011 ) have suggested that most 



of the ISO emission comes from non-dissociative C-shocks with shock velocities of order 15 
to 25 km s~^ and pre-shock densities of order 10^ to 10^ cm~^. A significant enhancement of 
water in these post-shocked regions is needed to explain the ISO observations. Many note 
that the strong [01] emission detected in all of these outflows and the d etection of [Sill] iri 



some of the outflows, may require the presence of dissociating J-shocks. 



Benedettini et al. 



(120021 ) compared the water emission observed by SWAS and ISO in four outflows and 
suggested that warm, approximately 1000 K, gas was responsible for the ISO emission and 
some of the SWAS emission, however they also concluded that SWAS was also sensitive to 
a cooler gas component not traced by ISO. 



We have used the p ublished line 



2000), NGC1333-SVS13 flMohnari et al 



Ceccarelli et al. 



20021) an d Orion KL ( 



note that 



1998 ), HH25mm 



Lerate et al. 



Froebrich et al 



uxes for L144 8-mm (INisini et al 



2000), L1157 ( 



Benedettini et al. 



Giannini et al. 



2000), Ser SMMl flLarsson et al 



1999 



200lh. L1689N 



20061 ) to perform a simple C-shock model analysis. We 
(I2OO2I ) reported two pointings of ISO toward the Ceph A HW2 
outflow; however, II2O emission was not detected, so we have not included this source in 
our analysis. Our modeling is aimed at estimating the shock solid angle responsible for 
producing the ISO emission and whether any signiflcant fraction of the SWAS water line 
flux originates in this shocked gas. In Table 4 we present the flux in the 179.5 /im line 
of ortho-Il20 observed by ISO for the six outflows with detected emission that were in 
common with our survey. In L1448-mm, NGC1333-SVS13, and LI 157 multiple positions 
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were observed by ISO and we have combined the fluxes for positions that are within the 



SWAS beam. In addition to CO and H^O, ISO detected re 



emission in these outflows (Ceccare 



Molinari et al. 



2000 



Giannini et al 



i et a. 



2001 



1998 



Nisini et al 



Larsson et al. 



ative 



2OO2I). 



20001 : 



y strong 63 /im [0 1] 



Benedettini et al 



2000 



Kaufman fc Neufeldl (11996h . 



To analyze the ISO emission, we use the C-shock models of 
We have integrated the line intensity contributions in the post-shock gas of these C-shock 
models to the point where the gas has cooled to a temperature of 50 K. At that point the 
gas is too cool to contribute much intensity to the observed ISO lines and the freeze out 
of water onto dust grains may become important, which is not included in this model. In 
addition, the postshock gas has become indistingishable from the gas previously modeled. 
We can use the high-J CO and H2O lines in conjunction with these C-shock model results 
to constrain the shock velocity and pre-shock gas density. 

The flux in the high-J CO lines and in the 179.5 fim transition of ortho-H20 are 
used to establish the correct C-shock model. The intensity ratio of the high-J CO lines, 
specifically the ratio of the J=20-19 transition to J=15-14 transition, is very sensitive to 
the temperature of the post-shock gas, and therefore provides a constraint on the shock 
velocity. The ratio of the CO J=15-14 line flux to 179.5 /im water line flux is sensitive to 
pre-shock density. Although these line ratios vary somewhat from outflow to outflow, all 
outflows require shock velocities > 20 km s~^ to produce the observed ratio of CO lines. 
For example the average CO line flux ratio of the J=20-19/J=15-14 transitions observed is 
0.45, and this ratio requires a shock velocity of 40 km s~^ for a pre-shock gas density of 10^ 
cm~^ or 20 km s~^ for a pre-shock gas density of 10^ cm~^. To reproduce the observed flux 
ratio of CO J=15-14 to 179.5 yum H2O, typically about unity, requires post-shock densities 
of > 10^ cm"'^. We have adopted a C-shock model with a shock velocity of 25 km s^^ and a 
pre-shock density of 10^ cm~^. The column of post-shock H2 gas in this C-shock model to 
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the point where the gas has cooled to 50 K is 6 x 10^^ cm~^ and the integrated intensity 
of the 179.5 fim H2O hne produced by this column of post-shock gas is 8.1 x lO^'^ ergs s^^ 
cm~^ sr~^. The H2O abundance relative to H2 in the post-shock gas is 4 x 10~^. 

Based on this adopted C-shock model we estimate the total solid angle of shock front 
needed to reproduce the observe line flux and consequently the total mass of shocked gas; 
these results are given in Table 4. With the exception of Orion KL, the angular area of 
these shocked regions is only ~10 square arcseconds, much smaller that the angular area 
of the outflow mapped in the low-J CO lines or the SWAS and ISO beamsizes. However 



for Orion KL, the angular size is about 1280 square arcseconds, comparable to the angu 



size o f the outflo w mapped in the 



19991 ). Recently, 



Cernicharo et al 



ow-J rotational lines of CO (iRodriguez-Franco et al. 



ar 



( I2OO6I ) modeled over 70 far-infrared rotational lines of 
water observed by ISO with the LWS. They modeled the outflow as a 40" diameter shell 
expanding at 25 km s^^ with a temperature between 80-100 K, density of 2.5- 3.5x10^ cm~ ^ 



and a relative abunda nce of water of 2- 3 x10 ^ , smaller than that derived by 



Harwit et al 



mm . The resuhs of 



Cernicharo et al. 



(120061 ) are similar to what we quoted earlier based 



on the SWAS data. Unlike the other outflows, the water emission observed by ISO in Orion 
KL may arise from the same gas that is responsible for producing the low-J CO emission. 

With exception of Orion KL, the mass of warm gas needed to reproduce the ISO 
observed lin es is only 0.0008 t o 0.005 solar masses. Our mass estimates are similar to those 



reported by 



Benedettini et al. 



fl2002l ) based on a thermal model with temperatures between 
600 and 1400 K. Thus, the mass of warm gas responsible for producing the ISO lines is only 
about 1% of the total outflowing gas mass inferred from our earlie r analysis of the J=l-0 



12 



CO line. For Orion KL, the mass in the shell model presented by 



Cernicharo et al. 



fl2006h 



is of order 6 solar masses, similar to our total mass estimate from the low-J CO lines. 



These warm, postshock gas regions probed by ISO also produce emission in the 538 fim 
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line observed by SWAS. The adopted C-shock model has a flux ratio of the 179 //m/538 /im 
lines of ortho-H20 of approximately 32. The observed flux ratio of these lines is summarized 
in Table 4. It is important to remember that the ISO observations are made in a smaller 
beam than SWAS (~75"), so that a single ISO observation may not detect all of the 179 
/im emission that is contained within a SWAS beam. Even in the outflows with multiple 
ISO pointings, the entire SWAS beam is not fully sampled, thus the observed line flux 
ratios in Table 4 are only lower limits. We estimate that the fraction of the SWAS water 
line flux that originates in th e warm, post-shock gas t o be as little as 10% in some outflows 



to more than 50 % in others. 



Benedettini et al. 



(j20o3) reached a similar conclusion. It may 



be signiflcant that the two outflows with the most complete ISO coverage (L1448-mm and 
Orion KL) have the largest line flux ratios. A signiflcant fraction of the SWAS water line 
flux may be arising from the very small mass of warm post-shock gas. 



4.7. Abundance of [CI] in Outflows 



SWAS obtained very high signal to noise spectra of [CI] toward these outflow sources. 
Obvious high velocity [CI] emission was detected in Mon R2, L1228, Ceph A IIW2, Orion 
KL, and S140 and spectra of [CI] in these sources are shown in Figure 2. Although no 
high velocity emission was detected in L1448-mm, the baseline in this source was somewhat 
worse. Also in Ser SMMl, the blue wing of the [CI] line was contaminated by emission in 



the reference position. T 
was previously noted by 



le presence o f high 



Walker et al. 



veloc ity [CI] emissio n in rn any of these outflows 



fll993h and 



Minchin et al. 



fll994h . 



In Orion KL, the spectrum shown in Figure 2 appears to have an additional emission 
feature at a velocity of about 45 km s~^. It is unlikely that this feature is related to the 
redshifted outflowing gas, as the line proflle would have to be signiflcantly different from 
those observed in either CO or H2O. The line is unlikely to arise from the opposite sideband. 
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Since SWAS does not doppler track in real time, frequency corrections are applied to each 
short time segment in the data pipeline, and thus lines in the opposite sideband should not 
add coherently to produce such a discrete feature. Assuming this feature is in the same 
sideband as [CI] and originates in gas with a Vlsr of 9 km s~^, then the line frequency is 
approximately 492.10 GHz. We searched the JPL line list for relatively low excitation lines 
of molecules known to be present in Orion and the only candidates are the 133^io-132,ii 
transition of CHsCHO-e at 492.093 GHz and the 138,6-137,6 transition of C^H sOH at 492.08 6 



GHz . The lines of CH.sCHO are ex tremely weak in the spectrum of Orion (jSutton et al 



1985 



Ziurys and McGonagle 



19931), and thus is probably an unlikely identification for the 



feature seen by SWAS. Ethanol, on the other hand has a number of lin es detected in the 



19951 ). but the frequency 



336 to 340 GHz range at several locations in Orion (iSutton et al.l 
of the transition is significantly lower than is needed to explain the SWAS feature. None of 
the previously published spectra of [CI] in Orion KL had adequate sensitivity to confirm or 
refute this feature. 

Using the same velocity intervals as defined earlier for the redshifted and blueshifted 
emission in each outfiow, we have computed the [CI] outfiow line fiux. A summary of the 
[CI] line fiux is presented in Table 5. We note that we have included the feature seen within 
the redshifted emission in Orion KL as part of the integrated intensity of the outfiow. 
Although we have [CI] detections in these velocity intervals for most of the outfiow sources, 
in many cases the emission may be dominated by the gaussian wings of the quiescent line 
emission. We computed the CI abundance in an identical manner as we computed the 
H2O abundance, using the same physical model for each of the outfiows. The resulting CI 
abundance relative to H2 is summarized in Table 5. 

We find that the abundance of CI relative to H2 is typically about 1x10^^. The largest 
abundance of CI is found for NGC2264 D and p Oph A. However, in NGC2264 D it is 
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questionable whether the weak [CI] emission in the red or blue shifted velocity intervals is 
related to the outflow. Very weak (T^ ~ 0.05 K) and very broad [CI] emission is detected 
toward p Oph A, similar to the emission profile seen in H2O. The CI abundance in p Oph 
A is about five times larger than in the other outflow. For quiescent cloud emission the 



aburidance of CI relativ e to CO is typically 0.1 to 0.5 (IZmuidzinas et al. 



2000 



1988 



Plume et al. 



Howe et al. 



2OOOI ). similar to that determined for the outflowing gas if we assume a 



CO to H2 ratio of 1x10 ^. Thus the CI in the outflowin g gas is simila r in ab undance to 



Walker et al. 



fll993h . The [CI] 



that of the ambient gas in agreement with the results of 
spectra have insufficient signal to noise to investigate variations of the CI abundance with 
outflow velocity. 



5. Discussion and Conclusions 

The lio — loi transition of ortho-H20 at 538 /xm has been detected by SWAS in 18 
molecular outflows. The II2O line proflles are similar to the line proflles observed for 
the J=l-0 transition of ^^CO and suggest that the emission seen in both species may be 
produced by the same gas. If we assume that the SWAS II2O emission arises in the same 
gas that makes up the bulk of the molecular outflow, we flnd that the outflowing gas 
has an ortho-H20 abundance relative to II2 typically between 10~® and 10^^. However, 
there are a few exceptions: most notably Orion KL and LI 157 have anomalously high 
relative ortho-H20 abundances of about 10~^, and Mon R2 has an anomalously low relative 
abundance of about 10~^. The relative ortho-H20 abundances in Table 2 have substantial 
uncertainty that arise almost solely due to the sensitivity of the derived water abundances 
to the assumed temperature and density of the outflowing gas, both of which are poorly 
known. The ^^CO J=5-4 observations by SWAS suggest that our physical model cannot 
be too much in error, and thus we believe that the abundance of water in most of these 
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outflows is elevated relative to that measured in quiescent cloud gas. We also derive the 
CI abundance in the outflowing gas, and flnd values that are similar to quiescent cloud 
material, and thus unlike water, the abundance of atomic carbon appears to be unaffected 
by the outflow activity. 

We also analyzed the velocity dependence of H2O, and find that the abundance of 
ortho-H20 varies significantly with velocity. In nearly all of the outfiows, we find a steady 
increase in the ortho-H20 abundance with increasing radial velocity of the outfiowing 
gas. In the most striking examples (NGC 1333 SVS13 and HH25mm), the abundance 
of water increases by two orders of magnitude from the lowest velocity to the highest 
velocity outfiowing gas. Nevertheless, the mass of outfiowing gas with greatly elevated H2O 
abundance is very small, representing at most only ~1% of the total outfiow gas mass. 

The ISO observations provide evidence for the presence of a much warmer outfiow 
component in those outfiows observed. With the exception of Orion KL, our C-shock 
modeling of the ISO emission suggest that this warm gas constitutes less than 1 % of the 
total outfiow mass and arises from an extremely small fraction of the molecular outfiow 
solid angle. However, this gas is expected to have a relative abundance of water in excess of 
10~^. The presence of this warmer component adds an additional complexity to abundance 
determinations in the cooler gas. Our shock modeling, although simple, suggests that a 
significant fraction of the water fine fiux observed by SWAS could arise in this warmer gas 
component. Unfortunately the ISO observations did not cover the full extent of the SWAS 
beam nor did it have velocity resolution, so it is extremely uncertain what fraction of the 
SWAS line fiux is associated with this warm gas probed by ISO. We estimate that for the 
outfiows with both ISO and SWAS observations, as much as 50 % of the SWAS fiux could 
originate in the warmer gas component. If any significant fraction of the SWAS H2O line 
fiux arises from the warm gas, then we have overestimated the line fiux from the cooler 
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gas traced by the low-J transitions of CO, and consequently overestimated the ortho-H20 
abundance in this cool gas. Therefore our estimates of the water abundance presented in 
Table 2 are likely to be too large. The extent to which the abundance is overestimated 
is very uncertain, but if only half of the SWAS line flux is due to the warm gas this only 
reduces the H2O abundance by a factor of two, which is well below other uncertainties in 
our abundance determination. Thus, unless nearly all of the SWAS line flux is due to this 
small mass of warm, shocked gas, the ortho-H20 abundance in the low velocity outflowing 



gas is enhanced relative to tha t measured 
only 10~^ to 10"*^ according to 



Snell et al. 



br qui escent cloud gas (relative abundance of 



f 2000bl )). 



The studies by 



Kaufman fc Neufeldl (Il996l ) and 



Bergin et al. 



(Il998l ) demonstrated that 



even a mild C-shock with shock velocities greater than 10 km s ^ will elevate the gas 
temperature to permit the efficient conversion of oxygen into water and to vaporize a large 



fraction of the water ice on grains. Thus, most of the oxygen, not in CO, wil 
converted into water resulting in relative water abundances in excess of 10~^. 



be rapidly 



Bergin et al. 



(119981 ) suggested that these elevated levels of water will persist long after the gas has cooled 
and any enhancement of the water abundance may be present over the entire lifetime of 
these outflows. The relative abundance of H2O is thus a good tracer of the dynamical 
history of the outflowing molecular gas. Our analysis of the SWAS data rules out the 
possibility that most of the outflow gas has passed through a shock in excess of 10 km s~^. 

The ISO and high velocity SWAS results for some outflows do indicate highly elevated 
water abundances in a small fraction of the outflowing gas mass that is consistent with 
moderate velocity C-shocks. The current observations however do not have sufficient spatial 
and/or velocity resolution to determine exactly how these emissions might be related. It 
is suggestive that since the mass associated with the ISO and high velocity SWAS II2O 
emission are very similar that these may arise from the same C-shock. For these outflows it 
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would be extremely interesting to determine the spatial and kinematic relation between the 
ISO emission, the high velocity SWAS emission and well-established tracers of C-shocks, 
such as the 2.12 /xm H2 emission seen in ma ny of these outflows such as L1448-mm 



(IDavis fc Smith 



I995I) and NGC1333-SVS13 flMolinari et al.ll2000h . Future observations 



with Herschel may help out in this regards. 



by 



Mechanisms for suppressing the H2O abundance in the postshock gas were explored 



Snell et al 



(l2005l ). The only mechanism that is probably relevant for these outflows is 



the freeze-out o 
in the model of 



gas-phase \v a ter on to grain mantles. This was the primary mechanism 



Bergin et al. 



(119981 ) for the reduction of gas-phase water in post-shock 



gas. The timescale for signiflcant depletion of gas-phase water depe nds on the po s t shoc k 



Bergin et al 



mm 



density of the gas and the total grain cross-section. In the model of 
for gas denser than 10^ cm~'^, the timescale is shorter than 10^ years, however at densities 
of 10^ cm~^ the time scale is longer than 3x10^ years. However since outflow lifetimes are 
extremely poorly understood, it is possible that the freeze-out o: 
factor. Based on the timescale for water freeze-out 



Snell et al. 



water is an important 
(120051 ) provided a simple 



expression for the hydrogen column density in the post-shock gas to the point where water 



will be frozen onto grain man tles. This ca 



than the numerical model of 



Bergin et al 



culation uses a slightly larger grain cross-section 



( 1l998l ). Assuming a shock velocity of 10 km 



s^^ and a shock compression factor of 10, the molecular hydrogen column would be of 
order 10^^ cm~^. If the relative abundance of ortho-Il20 in this column was as >10~^, as 
predicted by the C-shock models, the emission from this gas would greatly exceed what 
SWAS observed. Unless a more efficient mechanism is found for depleting the gas-phase 
water in the post-shock gas, we believe that the H2O abundance in outflows provides strong 
evidence that the bulk of the outflowing molecular gas was never subjected to shocks 
greater than 10 km s~^ that would convert all of gas phase oxygen to water. This result 
places severe restrictions on the mechanism by which the molecular gas is accelerated in 
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these outflows and favor meclianisms, sucfi as turbulent entrainment, tliat can accelerate 
the gas gently. However, some enhancement of the water abundance in the bulk of the 
outflowing gas over that measured in quiescent cloud gas may be necessary. Slower shocks 
(< 10 km s~^) could heat the gas and dust to temperatures above the water-ice evaporation 
temperature, but below the te mperature in which the r apid neutral-neutral reactions are 



(120061 ) for the water abundance measured 



activited, a suggestion made by iCernicharo et al. 
in Orion KL.. 

Finally we note that strong [01] emission at a wavelength of 63 fim was observed by 
ISO in all of the outflows. Many of the papers presenting ISO outflow results suggest that 
the origin of the [01] emission is in J-shocks that dissociate the molecular gas. Alternatively, 
the [01] emission could arise from the same weak shocks that accelerate the bulk of the 
molecular gas. Future observations with Herschel, which has better angular and spectral 
resolution, may help determine the relationship between the H2O and [01] emissions and 
other shock tracers in these outflows and provide a better understanding of the evolution of 
the H2O abundance in these outflows. 
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Table 1. Positions of Observed Outflows & SWAS Observing Times 



Source a(2000)i 5(2000)^ Distance Unt 

pc hr 



J_/i44o-mm 


OQ 


zo 


oU. 


.0 


1 on 
+oU 


A K 


A Q 


oUU 


Q 1 /I 7 

oi.4^ 


r^j A on 


(Jo 


Zl 


QQ 
OO, 


r 
.0 


+00 


Atl 

4o 


00 


yuu 


on rn 

oy.ou 


i\LiLyiooo-o Vbio 


Do 


on 


Uc), 


. 1 


1 Q1 

+ol 


1 

io 


(Jo 


ooo 
ZZU 


on 

zy.oz 


Orion KL 


Do 


00 


1 /I 

14, 


r 

.0 


-Uo 


OO 


O / 


oUU 


o.ov 


0MC2 


05 


35 


27, 


.3 


-05 


09 


49 


450 


6.34 


HH25nim 


05 


46 


07, 


,3 


-00 


13 


40 


400 


48.02 


NGC2071 


05 


47 


04, 


,1 


+00 


21 


43 


390 


19.32 


Men R2 


06 


07 


46, 


,7 


-06 


22 


42 


950 


17.55 


NGC2264 D 


06 


41 


03, 


.9 


+09 


34 


39 


800 


35.55 


NGC2264 C 


06 


41 


10. 


.7 


+09 


29 


07 


800 


11.36 


p Oph A 


16 


26 


23, 


,4 


-24 


23 


02 


160 


59.68 


L1689N 


16 


32 


22, 


,7 


-24 


28 


33 


120 


24.12 


Ser SMMl 


18 


29 


49, 


.6 


+01 


15 


20 


310 


35.60 


L1157 


20 


39 


06, 


,5 


+68 


02 


14 


440 


55.45 


L1228 


20 


57 


13, 


.0 


+77 


35 


47 


300 


33.87 


IC 1396N 


21 


40 


42, 


.3 


+58 


16 


10 


750 


59.02 


S140 


22 


19 


17, 


.1 


+63 


18 


46 


910 


24.78 


Cepheus A HW2 


22 


56 


17. 


.9 


+62 


01 


50 


725 


17.62 



^ Units of right ascension are hours, minutes, and seconds, and 
units of decUnation are degrees, arcminutes, and arcseconds. 
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Table 2. SWAS Line Fluxes and H2O Abundances^ 



Source 




Vlsr Interval 


fTX{H20)dv 


H2O Line Flux 


12QO/13CO 


Mass 


0-H2O Abundance'' 






(km s 




(K km s-l) 


(10-20 w cm-2) 




(Mm) 




L1448-mm 


Blue 


-31.0 


-1.0 


0.28 


(0.05) 


0.75 (0.13) 


21.4 ( 4.0) 


4.9x10-1 


1.5x10-'^ 




Keel 


8.0 


23.0 


0.27 


(0.03) 


n '70 /'n An\ 

U.7z (u.uyj 


21.6 ( 4.0) 


n^^in — 1 

z.yxiu 


O./XlU 


GL490 


Blue 


-42.0 


-17.0 


0.57 


(0.03) 


1.51 (0.09) 


19.6 ( 0.7) 


1.4x10^ 


1.3x10-*^ 




Keel 


-8.0 


22.0 


0.34 


(0.04) 


n on /n 1 r\\ 
U.yU (^U.lUj 


27.8 ( 1.7) 


1.2 X 10"*^ 


y.oxiu 


NGC1333-SVS13 


Blue 


-20.5 -» 


4.5 


0.54 


(0.04) 


1.42 (0.11) 


48.3 ( 9.3) 


5.4x10-1 


1.6x10-^ 




Red 


11.5 -» 


31.5 


0.64 


(0.04) 


1.08 (0.10) 


33.4 ( 3.1) 


6.7x 10— 1 


1.3X10 


Orion KL 


Blue 


-35.5 -» 


4.5 


40.5 


(0.09) 


106.43 (0.25) 


27.6 ( 0.8) 


1.2x101 


9.3x10-6 




Keel 


13.5 -» 


53.5 


48.1 


(0.09) 


1 Of; A^ { c\ oc^ 
lzo.41 


0*7 (\ t f\ n\ 

27. y ( u.yj 


9.8 


i.OXlU 


OMC-2 


Blue 


-2.5 


7.5 


0.84 


(0.07) 


2.21 (0.17) 


28.0 ( 3.5) 


3.4 


3.6x10-*^ 




Ked 


13.5 


23.5 




/n ^'7^ 
^U.U ( ) 


2.17 (U.17j 


62.1 (17.7) 


1.1 


0.7X1U 


HH25mm 


Blue 


-6.5 


8.5 


0.36 


(0.03) 


0.96 (0.09) 


21.0 ( 0.9) 


4.4 


4.2x10"'^ 




Red 


12.0 ^ 


37.0 


U.oU 




1.33 (0.12) 


25.2 ( 1.6) 


4.5 


7.0x10-"^ 




oiue 


-23.5 


6.5 


1.04 


(0.07) 


'7Q fC\ 1 '7^ 

z. i (^U.l * ) 


OK '7 ( C\ Q^ 

zo. / U.cS ) 


D.O 


1 1 \y 1 n — 7 
i.iXiU 




Red 


12.5 


37.5 


1.73 


(0.06) 


4.55 (0.16) 


33.4 ( 1.2) 


6.9 


2.1X10-''' 


MonR2 


rJlue 


-2.2 -+ 


7.8 


0.47 


(0.04) 


l.zo (O.iOj 


iO.3 ( U.lj 


1 Owl n2 
l.zXiO 


1.8XiO 




Red 


13.0 ->■ 


28.0 


0.59 


(0.05) 


1.56 (0.12) 


8.6 ( 0.1) 


2.6x10^ 


1.3x10-^ 


NGC2264 D 


Blue 


-9.0 ^ 


1.0 


0.13 


(0.03) 


0.33 (0.07) 




1.8 


1.4x10-"^ 




Red 


10.0 ->■ 


20.0 


0.19 


(0.03) 


0.51 (0.07) 




1.9 


2.4x10-6 


NGC2264 C 


Red 


11.0 


21.0 


0.83 


(0.04) 


2.18 (0.11) 


21.1 ( 1.1) 


7.2 


1.7x10-6 


p Oph A 


Blue 


-3.5 -+ 


1.5 


0.28 


(0.01) 


0.74 (0.04) 


17.8 ( 2.7) 


2.2x10-1 


9.4x10-'^ 




Red 


6.0 -> 


11.0 


0.17 


(0.01) 


0.43 (0.04) 


13.3 ( 1.1) 


4.3x10-1 


3.7x10-'^ 


L1689N 


Blue 


-12.5 


2.5 


0.57 


(0.05) 


1.49 (0.12) 


18.8 ( 2.4) 


2.7x10-1 


9.4x10-'^ 




Red 


6.0 -> 


21.0 


1.08 


(0.05) 


2.84 (0.12) 


20.8 ( 3.0) 


1.1x10-1 


2.9x10-6 
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Table 2 — Continued 



Source 




Vlsr Interval 


/T*(H20)di; 


H2O Line Flux 


12CO/ 




Mass 


0-H2O Abundance'' 






(km s 




(K km s-i) 


(10-20 W cm-2) 






(M©) 




Ser SMMl 


Blue 


-4.5 


5.5 


0.40 


(0.02) 


1.05 


(0.06) 


19.8 ( 


1.3) 


1.8 


7.1x10-'^ 




Red 


11.0 


21.0 


0.32 


(0.02) 


0.83 


(0.06) 


20.0 ( 


1.1) 


2.9 


3.8x10-'^ 


L1157 


Blue 


-13.5 ->■ 


1.5 


0.70 


(0.02) 


1.84 


(0.06) 


28.7 ( 


5.1) 


4.1x10-1 


8.0x10-'^ 




Red 


4.0 ->■ 


29.0 


0.60 


(0.03) 


1.58 


(0.08) 






3.2x10-1 


9.7x10-6 


L1228 


Blue 


-19.5 


-9.5 


0.17 


(0.02) 


0.45 


(0.06) 


21.7 ( 


1.3) 


8.7x10-1 


5.9x10-'^ 




Red 


-6.0 -> 


4.0 


0.14 


(0.02) 


0.38 


(0.06) 


20.7 ( 


2.5) 


5.1x10-1 


8.5x10-'^ 


IC1396N 


Blue 


-18.0 -+ 


-3.0 


0.36 


(0.02) 


0.95 


(0.06) 


31.3 ( 


6.0) 


1.5 


2.7x10-6 




Red 


4.0 -» 


19.0 


0.33 


(0.02) 


0.88 


(0.06) 


36.2 ( 


9.1) 


1.2 


3.2x10-6 


S140 


Blue 


-36.0 ^ 


-11.0 


0.24 


(0.04) 


0.63 


(0.11) 


15.7 ( 


0.4) 


3.2x101 


2.4x10-8 




Red 


-4.0 


6.0 


0.47 


(0.03) 


1.25 


(0.07) 


18.3 ( 


0.6) 


1.6x101 


1.0x10-''' 


Ceph A HW2 


Blue 


-35.0 


-15.0 


1.34 


(0.06) 


3.53 


(0.16) 


14.8 ( 


0.5) 


l.SxlQi 


1.9x10-'^ 




Red 


-6.0 


14.0 


1.27 


(0.06) 


3.35 


(0.16) 


19.0 ( 


0.9) 


1.7x101 


2.0x10-^ 



*The numbers in parantheses are l-cr statistical uncertainties on the integrated intensity, line flux and observed isotopic ratio 
l'o-H20 abundances relative to H2 assuming a H2/CO ratio of 10'* 
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Table 3. SWAS Line Fluxes and H2O Abundances as a Function of Velocity^ 



Source 


Vlsr Interval 


fT* (B.20)dv 


H2O Line Flux 




Mass 


0-H2O Abundance'^ 




(kms 


-1") 


(K km s-l) 


(10-20 w cm-2) 




(M©) 






-6.0 


-1.0 


\)X)i \\).\JZ) 


U.lo i^U.Uoj 


ol.o ty.o) 


l.D X lU 


1 9 ^ 1 n ~ 6 

l.Z X lU 




-11.0 


-6.0 


^ n Of! 
u.uo 


u.io 




9 1 V 1 
z. 1 X lU 


<^ 1 .1 X. lU 




-16.0 -» 


-11.0 


U.UD ^^U.Uz_) 


U.lo ^^U.Uo_) 




Q 'Tn/ 1 n— 3 
y. ( X lu 


1 9 \/ 1 n— 5 
l.z X lU 




-21.0 -» 


-16.0 


U.UD 


U.lo ^^U.UO_) 




Q /I \/ 1 n— 3 
y.4X lU 


l.oX lU 




-26.0 


-21.0 


<. U.UO 


< U.lo 




9 A \/ 1 n— 3 
z.OX lU 


^A c; \/ 1 n— 5 
<.4.o X lU 




-31.0 


-26.0 


^ n OR 

U.UD 


<, U.iD 




Q 9 V 1 n~3 

O.Z X lU 


<^^.^ X lU 




8.0 


13.0 


U. Iz ) 


U.ol (^U.Uoj 


oD.o y y.Uj 


1 Q ^ 1 n— 1 
1.0 X lU 


0.0 X lU 




13.0 


18.0 


U.Uo l^U.UZj 


U.Zl i^u.uoj 


y.u z."*^ 


0.0 X lU 


41. X lU 




18.0 -^ 


23.0 


n HQ i C\ AO\ 

U.UcS t^U.U^j 


n on ( c\ nc;\ 
U.zU l^U.Uoj 




A Q 1 n — 3 
4.0X lU 


1.1 X lU 


GL490 Blue 


-22.0 


-17.0 


0.15 (0.01) 


0.39 (0.04) 


19.1 ( 0.5) 


9.9 


4.8x10-^ 




-27.0 ^ 


-22.0 


0.11 (0.01) 


0.29 (0.04) 


21.2 ( 1.5) 


3.1 


1.1x10-*^ 




-32.0 


-27.0 


U.lo (^U.Ul^ 


U.o4 (^U.U4^ 


9C 9 f a '7\ 
ZO.Z 0. i J 


0. 1 X lU 


0.4X lU 




-37.0 -» 


-32.0 


U.UC5 ^^U.Ul ) 






1.0 X lU 


1.0 X lU 




-42.0 -» 


-37.0 


U.IU ^^U.Ul^J 






0.^/\ lU 


C5.0 X lU 


rtea 


-8.0 -» 


-3.0 


U.IU \^U.U1_) 


U.zo 


OA n f n Q\ 
z4:.u u.y_) 


0.0 


o.y X lu 




-3.0 ^ 


2.0 


U.Uo (^U.Ui ) 


n 00 ( c\ r\A\ 

V.ZZ [V.K)^) 


A A fi ^ V 0\ 
44. D [ i.Z) 


1 A 
1.4 


1. / X lU 




2.0 


7.0 


U.UO l^U.Ul^ 


U.l^ l^u.u^^ 


o't.o 1^1 1.0^ 


0.0 X lU 


9 7 V 1 
Z. I X lU 




7.0 — » 


12.0 


<. U.U41: 


--^ n 1 9 
< U.lz 




1 fi \/ 1 n— 1 

l.D X lU 


*C / .0 X lu 




12.0 -» 


17.0 


U.Ui ^^U.Ui^ 


U.lo \^U.U4l:^ 




o.oX lU 


9 \/ 1 n— 6 
z.o X lU 




17.0 -» 


22.0 


< 0.04 


< 0.12 




2.6xl0~2 


<4.5xl0-^ 


NGC1333-SVS13 Blue 


-0.5 ^ 


4.5 


0.20 (0.02) 


0.53 (0.05) 


34.9 ( 3.1) 


6.7x10-1 


5.1x10-''' 




-5.5 ^ 


-0.5 


0.13 (0.02) 


0.34 (0.05) 




2.6x10-2 


5.7x10-6 




-10.5 


-5.5 


0.09 (0.02) 


0.24 (0.05) 




5.5x10-3 


1.6x10-^ 




-15.5 -» 


-10.5 


0.12 (0.02) 


0.30 (0.05) 




2.4x10-3 


4.8x10-^ 




-20.5 ^ 


-15.5 


< 0.05 


< 0.14 




1.5x10-3 


<3.1xl0-^ 


Red 


11.5 ^ 


16.5 


0.30 (0.02) 


0.79 (0.05) 


30.3 ( 1.7) 


5.7x10-1 


7.2x10-'^ 




16.5 


21.5 


0.17 (0.02) 


0.44 (0.05) 




4.4x10-2 


4.6x10"^ 




21.5 ^ 


26.5 


0.10 (0.02) 


0.27 (0.05) 




8.7xl0~^ 


1.5x10-^ 




26.5 -» 


31.5 


0.07 (0.02) 


0.18 (0.05) 




1.2x10-3 


7.9x10-5 
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Table 3 — Continued 



Source Vlsh Interval /T*(H20)dj; H2O Line Flux i^^o/isco Mass 0-H2O Abundance'' 







(km s 




(K km s-i) 


(10-20 w cm-2) 




(Mo) 






Orion KL 


Blue 


-0.5 




4.5 


10.2 


(0.03) 


26.92 (0.09) 


20.2 ( 0.4) 


7.8 




2.4x10"'' 






-5.5 




-0.5 


8.28 


(0.03) 


21.77 (0.09) 


25.1 ( 1.1) 


2.2 




6.4x10-6 






-10.5 




-5.5 


6.30 


(0.03) 


16.57 (0.09) 


34.4 ( 3.1) 


1.2 




8.1x10-8 






-15.5 




-10.5 


4.98 


(0.03) 


13.11 (0.09) 


49.2 ( 9.5) 


5.3x10" 


1 


1.3x10-5 






-20.5 




-15.5 


3.80 


(0.03) 


9.99 (0.09) 




2.6x10" 


-1 


1.9x10-5 






-25.5 




-20.5 


3.01 


(0.03) 


7.92 (0.09) 




1.1x10" 


1 


3.6x10-5 






-30.5 




-25.5 


2.26 


(0.03) 


5.94 (0.09) 




3.8x10" 


2 


9.2x10^5 






-35.5 




-30.5 


1.60 


(0.03) 


4.21 (0.09) 




1.8x10- 


2 


1.2x10-"' 




Red 


13.5 


- 


18.5 


13.7 


(0.03) 


36.12 (0.09) 


24.1 ( 0.5) 


6.7 




4.3x10"® 






18.5 


~* 


23.5 


10.8 


(0.03) 


28.36 (0.09) 


27.2 ( 1.3) 


1.8 




1.3x10-5 






23.5 


~* 


28.5 


8.22 


(0.03) 


21.61 (0.09) 


28.3 ( 2.4) 


8.9x10" 


-1 


1.8x10-5 






28.5 


~* 


33.5 


5.94 


(0.03) 


15.61 (0.09) 




2.2x10" 


1 


4.9x10-5 






33.5 


~^ 


38.5 


3.99 


(0.03) 


10.50 (0.09) 




1.1x10" 


1 


6.6x10-5 






38.5 


~^ 


43.5 


2.61 


(0.03) 


6.87 (0.09) 




4.5x10" 


2 


9.6x10-5 






43.5 


~* 


48.5 


1.70 


(0.03) 


4.47 (0.09) 




2.8x10" 


2 


8.5x10-5 






48.5 




53.5 


1.09 


(0.03) 


2.86 (0.09) 




2.2x10" 


2 


8.1x10-5 


OMC2 


Blue 


2.5 




7.5 


0.46 


(0.05) 


1.22 (0.12) 


21.7 ( 1.6) 


4.2 




i.exio-'' 






-2.5 




2.5 


0.38 


(0.05) 


0.99 (0.12) 




1.9x10" 


2 


1.7x10-5 




Red 


13.5 




18.5 


0.56 


(0.05) 


1.46 (0.12) 


52.4 ( 9.7) 


1.3 




5.3x10-'^ 






18.5 




23.5 


0.27 


(0.05) 


0.71 (0.12) 




1.6x10" 


2 


1.3x10-5 


HH25mm 


Blue 


3.5 




8.5 


0.16 


(0.02) 


0.42 (0.05) 


17.8 ( 0.4) 


4.8 




1.7x10-'' 






-1.5 




3.5 


0.12 


(0.02) 


0.33 (0.05) 




7.6x10" 


2 


7.0x10-6 






-6.5 




-1.5 


0.08 


(0.02) 


0.21 (0.05) 




7.8x10" 


-3 


3.0x10-5 




Red 


12.0 




17.0 


0.19 


(0.02) 


0.49 (0.05) 


26.4 ( 1.1) 


3.6 




3.0x10-'' 






17.0 




22.0 


0.16 


(0.02) 


0.42 (0.05) 


37.3 ( 8.2) 


3.1x10" 


1 


3.8x10-'^ 






22.0 




27.0 


0.09 


(0.02) 


0.23 (0.05) 




3.4x10" 


2 


2.7x10-5 






27.0 




32.0 


< 0.06 


< 0.16 




8.2x10" 


3 


<1.0xl0-'* 






32.0 




37.0 


< 0.06 


< 0.16 




1.3x10" 


-3 


<7.2xl0-* 


NGC2071 


Blue 


1.5 




6.5 


0.35 (0.03) 


0.91 (0.07) 


20.9 ( 0.4) 


6.4 




4.2x10-* 



Table 3 — Continued 



Source Vtsii Interval /r*(H20)di) Ha O Line Flux i^qq/ISCO Mass 0-H2O Abundance^ 

(km s-l) (K km s"!) (IQ-^O W cm'^) (M©) 







-3.5 




1.5 


0.29 


(0.03) 


0.75 


(0.07) 


31.7 ( 2.0) 


8.9x10-1 


2.1x10-'^ 






-8.5 




-3.5 


0.17 


(0.03) 


0.43 


(0.07) 




9.0x10-2 


1.1x10-6 






-13.5 




-8.5 


< 


3.08 


< 


3.21 




2.6x10-2 


<1. 7x10-6 






-18.5 


-* 


-13.5 


0.16 


(0.03) 


0.42 


(0.07) 




1.6x10-2 


5.5x10-6 






-23.5 




-18.5 


< 


D.08 


< 


3.21 




Q 1 n — 3 
o.zXiU 


<1.4XiU 




Red 


12.5 




17.5 


0.81 


(0.03) 


2.14 


(0.07) 


2b. 4: [ U.4} 


7.2 


9.9x10 






17.0 




22.5 


0.41 


^0.03) 


1.09 


^0.07) 


cr\ 1 / n 1 \ 

bu.i ( y.ij 


A '7 1 n — 1 
4. / X iU 


6.1 x 10— 






22.5 




27.5 


0.25 


(0.03) 


0.65 


(0.07) 




8.8x10-2 


1.8x10-6 






27.5 




32.5 


0.16 


(0.03) 


0.42 


(0.07) 




3.2x10-2 


3.0x10-6 






32.5 




37.5 


0.10 


(0.03) 


0.26 


(0.07) 




8.5x10-3 


6.6x10-6 


MonR2 


Blue 


2.8 




7.8 


0.38 


(0.03) 


1.01 


(0.07) 


9.9 ( 0.1) 


1.3xl02 


1.5x10-8 






-2.2 




2.8 


0.08 


(0.03) 


0.22 


(0.07) 




1.4x10-1 


2.7x10-6 




Red 


13.0 




18.0 


0.45 


(0.03) 


1.19 


(0.07) 


8.0 ( 0.0) 


Z.4X10 


1.1x10—8 






18.0 





23.0 


0.09 


(0.03) 


0.24 


(0.07) 


17.9 ( 1.7) 


1.3x10^ 


3.1x10-8 






23.0 




28.0 


< 


3.08 


< 


3.21 




2.4x10-1 


<1.2xl0-6 


NGC2264 D 


Blue 


-4.0 




1.0 


0.09 


(0.02) 


0.24 


(0.05) 




1.6 


1.1x10-6 






-9.0 




-4.0 


< 


3.05 


< 


3.14 




5.2x10-2 


<1.6xl0-5 




Red 


10.0 




15.0 


0.13 


(0.02) 


0.35 


(0.05) 




1.8 


1.6x10-6 






15.0 




20.0 


0.06 


(0.02) 


0.16 


(0.05) 




7.5x10-2 


2.9x10-5 


NGC2264 C 


Red 


11.0 




16.0 


0.52 


(0.03) 


1.37 


(0.08) 


22.1 ( 1.0) 


6.3 


1.2x10-6 






16.0 




21.0 


0.31 


(0.03) 


0.80 


(0.08) 


14.9 ( 4.1) 


5.8x10-1 


6.4x10-6 


p Oph A 


Blue 


-3.5 




1.5 


0.28 


(0.01) 


0.74 


(0.04) 


17.8 ( 2.7) 


2.2x10-1 


9.5x10-'^ 




Red 


6.0 




11.0 


0.17 


(0.01) 


0.43 


(0.04) 


13.3 ( 1.1) 


4.3x10-1 


3.7x10-'^ 


L1689N 


Blue 


-2.5 




2.5 


0.37 


(0.03) 


0.97 


(0.07) 


18.3 ( 1.8) 


2.2x10-1 


7.5x10-'^ 






-7.5 




-2.5 


0.13 


(0.03) 


0.33 


(0.07) 




6.8x10-3 


8.4x10-6 






-12.5 




-7.5 


< 


3.08 


< 


3.21 




9.9x10-4 


<3.3xl0-5 




Red 


6.0 




11.0 


0.72 


(0.03) 


1.89 


(0.07) 


21.3 ( 2.5) 


9.0x10-2 


2.5x10-6 






11.0 




16.0 


0.25 


(0.03) 


0.66 


(0.07) 




4.7x10-3 


1.6x10-5 



Table 3 — Continued 



Source Visii Interval /r*(H20)di) Ha O Line Flux i^qq/ISCO Mass 0-H2O Abundance'' 







(km s 




(K km s-i) 


(10-20 W cm-2) 




(Me) 








16.0 




01 n 


nil /a r\Q\ 
U.il (U.UoJ 


n on ( c\ f\'7\ 
U.zy l^U.U / ) 




A n 1 A — 4 
4.y X iU 


Q 1 1 n 
cS.i X iU 


Ser SMMl 


Blue 


0.5 


-* 


5.5 


0.27 (0.01) 


0.70 (0.04) 


17.8 ( 0.8) 


1.9 


4.7x10 






-4.5 




0.5 


0.13 (0.01) 


0.35 (0.04) 




1.9x10-2 


1.6x10 




rteci 


11.0 




1 n 
iO.U 


A 1 c /A ni A 
U. lo (U.Ui ) 


A /I Q (C\ (\A\ 

U.4c5 l^U.U4j 


1 7 n f' n fl^ 

i / .y u.D; 


Q A 
O.U 


1 1 n 
z.l X iU 






16.0 




01 n 


A "1/1 /A^1^ 
U.14 (U.Ui ) 


n Qft /n A/i\ 
U.oO (^U.U4j 




Q ^ 1 A — 2 
O.DX iU 


1 A r\ 
i.4X iU 


L1157 


Blue 


-3.5 




1.5 


0.33 (0.01) 


0.86 (0.04) 


35.4 ( 5.6) 


2.8x10-1 


5.3x10 






-8.5 




-3.5 


0.22 (0.01) 


0.58 (0.04) 




1.9x10-2 


4.9x10 






-13.5 




-0.0 


A 1 c; /a ni ^ 
U.iO (U.Ul ) 


n /in f'n (\a\ 
U.4U l^U.U4j 




Q n 1 A — 3 
o.y X iU 


Q c w 1 n 
0.0 X iU 




rtea 


4.0 




y.u 


A OQ /^A A1 ^ 
U.Zo \U.\Ji-) 


n c;n /n A/1^ 
U.oy (^U.U4^ 


QQ 7 / 7 A^ 

00. ( \ (.KJ) 


1 A — 1 

z.y X iu 


Q Q 1 n 
o.oX lU 






9.0 




14. U 


r\ 1 A (C\ A1 ^ 
U.i4 (^U.UiJ 


n /n A/1^ 
U.o* (^U.U4^ 




c; \^ 1 A — 2 
O.ZX iU 


1 Q 1 n 
i.oX iU 






14.0 




inn 


All (c\ ni ^ 
U.ii ^^U.UiJ 


n OQ /n A/1 ^ 
U.zcS 




Q R 1 A — 2 

O.DX iU 


1 c; 1 n 
i.O X lU 






19.0 





OA n 


A Afi ((\ ni ^ 
U.UD l^U.Ul ) 


C\ Af< (C\ A/1 ^ 
U.iD l^U.U4j 




n V 1 A— 2 

z.U X iU 


1 fi \/ 1 n 
i.D X iU 






24.0 




29.0 


0.07 (0.01) 


0.18 (0.04) 




7.0x10-3 


5.2x10 


T 1 OOQ 


rSiue 


-14.5 




-y.o 


A nn /"a ni ^ 
U.Uy ^^U.Uij 


U.z4 (^U.U4j 


OA n / 1 A^ 

zU.y i.Uj 


Q A A — 1 


Q 1 n 
O.ZX lU 






-19.5 




1 /I c; 
-14. 


A AQ / A ni ^ 

U.Uo (U.Ui ) 


A 00 ( C\ (\A\ 
V.ZZ [V.K)^) 




\y 1 A — 2 
z.Z X iU 


1 1 n 
i.z X iU 




-TLea 


-6.0 




-i.U 


A AC / A A1 ^ 

U.Uo (U.Ui ) 


A 01 (r\ A/i 'i 
V.Zl (^U.U4j 


01 1/0 A^ 

Zl.L Z.U ) 


c; 1 \y 1 A— 1 
O.i X iU 


A c; \^ 1 n 
4.0 X iU 






-1.0 






A Cil (c\ ni ^ 

U.U/ l^U.Ul J 


n 1 7 /n A/1 ^ 
U.l( l^U.U4j 




C; Q \y 1 A — 3 
u.oX iU 


A A \^ in 

4.4 X lU 


IC1396N 


Blue 


-8.0 




-3.0 


0.17 (0.01) 


0.46 (0.04) 


30.9 ( 4.4) 


1.2 


1.6x10 






-13.0 




-8.0 


0.10 (0.01) 


0.27 (0.04) 




1.1x10-1 


1.2x10 






-18.0 




-13.0 


0.08 (0.01) 


0.22 (0.04) 




1.6x10-2 


5.8x10 




Red 


4.0 




9.0 


0.20 (0.01) 


0.52 (0.04) 


30.3 ( 4.4) 


1.4 


1.7x10 






9.0 




14.0 


0.10 (0.01) 


0.26 (0.04) 




3.6x10-2 


2.6x10 






14.0 




19.0 


< 0.04 


< 0.11 




4.5x10-3 


<8.5xl0 


S140 


Blue 


-16.0 




-11.0 


0.14 (0.02) 


0.37 (0.05) 


13.8 ( 0.2) 


3.0x10^ 


1.5x10 






-21.0 




-16.0 


< 0.05 


< 0.14 


21.0 ( 1.3) 


3.8 


<4.1xlO 






-26.0 




-21.0 


< 0.05 


< 0.14 




1.6x10-1 


<8.7xl0 






-31.0 




-26.0 


0.09 (0.02) 


0.25 (0.05) 




2.8x10-2 


9.0x10 






-36.0 




-31.0 


< 0.05 


< 0.14 




1.5x10-2 


<1.5xl0 
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Table 3 — Continued 



Source 




Vlsr Interval 


/T*(H2 0)d« 


H2O Line Flux 


12QO/13CO 


Mass 


0-H2O Abundance'' 






(km s' 




(K km s-l) 


(10-20 W cm-2) 




(Mo) 






Red 


-4.0 


1.0 


0.42 (0.02) 


1.09 (0.05) 


17.3 ( 0.5) 


1.5x10^ 


9.7x10-^ 






1.0 -+ 


6.0 


0.06 (0.02) 


0.15 (0.05) 




4.7x10-1 


4.9x10-'^ 


Ceph A HW2 


Blue 


-20.0 


-15.0 


0.62 (0.03) 


1.63 (0.08) 


12.0 ( 0.2) 


1.9x10^ 


8.5x10-* 






-25.0 -+ 


-20.0 


0.34 (0.03) 


0.90 (0.08) 


29.0 ( 5.3) 


8.8x10-1 


8.5x10-'^ 






-30.0 -> 


-25.0 


0.25 (0.03) 


0.65 (0.08) 




1.1x10-1 


4.7x10-6 






-35.0 -> 


-30.0 


0.13 (0.03) 


0.35 (0.08) 




4.6x10-2 


5.5x10-6 




Red 


-6.0 


-1.0 


0.74 (0.03) 


1.94 (0.08) 


15.3 ( 0.4) 


1.9x10^ 


1.1x10-'^ 






-1.0 


4.0 


0.40 (0.03) 


1.04 (0.08) 




4.4x10-1 


2.1x10-6 






4.0 


9.0 


< 0.09 


< 0.24 




3.1x10-1 


<8.6xl0-'^ 






9.0 


14.0 


< 0.09 


< 0.24 




2.8x10-2 


<6.4xl0-6 



*The numbers in parantheses are 1-cr statistical uncertainties on the integrated intensity, line flux and observed isotopic ratio. All 
limits quoted are 3-it 

l'o-H20 abundances relative to H2 assuming a H2/CO ratio of 10* 
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Table 4. ISO Results 



Source 


Flux(o-H20 179 fj,m) 


Solid Angle 


Mass 


H2O Flux Ratio (179 /ini/538 nm) 




10-20 W cm-2 


sr 


Mq 




L1448-mm 


26.9±3.4^ 


4.3x10-1° 


4x10-3 


18.3 


NGC1333-SVS13 


14.7±4.lb 


1.8x10-1° 


8x10-^ 


4.7 


HH25mm 


10±2 


1.2x10-1° 


2x10-3 


4.4 


L1689N 


21±2 


2.6x10-1° 


3x10-^ 


4.8 


Ser SMMl 


14.0±3.2 


1.7x10^1° 


2x10-3 


7.4 


L1157 


26.5±3.5'= 


2.7x10-10 


5x10-3 


7.7 


Orion KL 


2550±300 


3.1x10-* 


5x10-1 


10.7 



^Includes positions L1448-mm and IRS3 
''Includes positions SVS13 and HH7 
■^Includes positions ON and Blue 
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Table 5. SWAS CI Abundances^ 



Source 



jTX{[Cl])dv CI Abundance'' 



(K ku 



L1448-mm 



GL490 



Blue 
Red 



< 0.14 

< 0.10 



Blue 0.49 (0.06) 
Red 0.53 (0.06) 



NGC1333-SVS13 Blue 0.55 (0.05) 
Red < 0.14 



Orion KL 

OMC-2 

HH25mm 

NGC2071 

MonR2 

NGC2264 D 

NGC2264 C 
p Oph A 

L1689N 



Blue 1.16 (0.03) 
Red 2.04 (0.03) 



Blue 
Red 



Red 



< 0.21 

< 0.21 



Blue 0.44 (0.04) 

Red 0.49 (0.05) 

Blue 0.49 (0.06) 

Red 0.49 (0.06) 

Blue 2.27 (0.05) 

Red 3.80 (0.07) 

Blue 0.26 (0.03) 

Red 0.29 (0.03) 



< 0.20 



Blue 0.80 (0.02) 

Red 0.83 (0.02) 

Blue 0.27 (0.05) 
Red < 0.15 



<1. 5x10-5 
<2.5xl0-5 

1.3x10-5 
1.7x10-5 

3.2x10-5 
<5.7xl0-s 

6.7x10-8 
1.4x10-5 

<9.7xl0-6 
<2.4xl0-5 

1.0x10-5 
1.3x10-5 

6.0x10-6 
6.5x10-6 

1.0x10-5 
8.9x10-6 

5.9x10-5 
7.2x10-5 

<8.5xl0-6 

5.5x10-5 
3.6x10-5 

9.0x10-6 
<8.1xl0-6 
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Table 5 — Continued 



Source 




/Tl([CI])d« 
(K km s-l) 


CI Abundance'' 


Ser SMI 


Blue 


< 0.12 


<4.2xl0" 


-6 




Red 


0.39 (0.04) 


9.0x10" 


-6 


L1157 


Blue 


0.10 (0.03) 


2.3x10" 


-5 




Red 


0.11 (0.04) 


3.4x10" 


-5 


L1228 


Blue 


0.31 (0.02) 


2.1x10" 


-6 




Red 


0.28 (0.02) 


3.3x10" 


-6 


IC1396N 


Blue 


< 0.13 


<2.0xl0" 


-6 




Red 


< 0.13 


<2.6xl0" 


-6 


S140 


Blue 


0.88 (0.03) 


1.0x10" 


-6 




Red 


0.37 (0.02) 


9.2x10" 


-6 


Ceph A HW2 


Blue 


1.36 (0.04) 


2.2x10" 


-5 




Red 


l.^ie (0.04) 


2.3x10" 


-5 



*The numbers in parantheses are l-a statistical uncer- 
tainties on the integrated intensity, line flux and observed 
isotopic ratio. All limits quoted are 3-cr 

''CI abundances relative to H2 assuming a H2/CO ratio 
of 10* 




Fig. 1.— (a) Spectra towards L1448-mm, GL490, NGC1333-SVS13, OMC 1, OMC 2 and 
HH25-mm. The positions observed are provided in Table 1. For each source there are 
two nanels. the imner nane] shows the snectrum of the Im — Im transition of ortho-HoO 




Fig. 1.— (b) The same as Fig. la for NGC2071, MonR2, NGC2264 D, NGC2264 C, p Oph 
A, and L1689 N. 
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Fig. 1.— (c) The same as Fig. la for Ser SMMl, L1157, L1228, IC 1396N, S140 and Ceph 
A HW2. 
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Fig. 2.— Spectra of the ^Pi ^Pq transition of [CI] towards OMC 1, Mon R2, L1228, S140 
and Ceph A HW2 obtained with SWAS. The positions observed are provided in Table 1. 



